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We report on a study of intrinsic superconductivity in a Weyl metal, i.e. a doped Weyl semimetal. 
Two distinct superconducting states are possible in this system in principle: a zero-momentum 
pairing BCS state, with point nodes in the gap function; and a hnite-momentum FFLO-like state, 
with a full nodeless gap. We hnd that, in an inversion-symmetric Weyl metal the odd-parity BCS 
state has a lower energy than the FFLO state, despite the nodes in the gap. The FFLO state, on 
the other hand, may have a lower energy in a noncentrosymmetric Weyl metal, in which Weyl nodes 
of opposite chirality have different energy. However, realizing the FFLO state is in general very 
difficult since the paired states are not related by any exact symmetry, which precludes a weak- 
coupling superconducting instability. We also discuss some of the physical properties of the nodal 
BCS state, in particular Majorana and Fermi arc surface states. 


I. INTRODUCTION 

The study of the interplay of superconductivity and 
nontrivial electronic structure topology has a long his¬ 
tory, dating back to the work on superfluid The 

recent discovery of topological insulators (Tlj^has rein¬ 
vigorated this subject. Proximity-induced superconduc¬ 
tivity on a 3D TI surface has been proposed as a route 
to realize Majorana fermions,^ and bulk topological su¬ 
perconductivity has also been studied extensively.^^^f^ 

Most recently, Weyl, and closely related Dirac 
semimetals, first discovered theoreticalljA^^^ and lately 
realized experimentally, E^HUI have come to the forefront 
of research on topologically-nontrivial phases of mat¬ 
ter. Perhaps the most interesting new feature that Weyl 
semimetals bring to the subject is that they are gap¬ 
less. The realization that gapless states of matter may 
be topologically nontrivial, just as gapped insulators, is 
of signihcant importance. 

The nontrivial electronic structure topology of Weyl 
semimetals arises from points of contact of nondegener¬ 
ate conduction and valence bands, which act as monopole 
sources of Berry curvature and thus carry an integer 
topological charge. A semimetal is realized when the 
Fermi level coincides with the points of contact and the 
valence band is filled while the conduction band is empty. 
This situation, however, is nongeneric and is unstable to 
impurities which will always give rise to unintentional 
doping, shifting the Fermi level into the conduction or 
valence bands. This motivates the study of Weyl metals, 
i.e. lightly-doped Weyl semimetals. It turns out that 
Weyl metals share most of the topologically-nontrivial 
properties with undoped Weyl semimetals.l^S'^Two fea¬ 
tures of the Weyl metals, which are directly related to the 
Weyl node topology, are very important in this regard. 
One is that, when the Fermi energy is close enough to 
the Weyl nodes, the Fermi surface breaks up into discon¬ 
nected sheets, each enclosing one Weyl node. The flux 
of the Berry curvature through each such Fermi surface 
sheet is equal to the topological charge of the correspond¬ 


ing node, which endows each 2D Fermi surface sheet with 
a nontrivial topological invariant, a Chern number. The 
second important property of a Weyl metal is the linear¬ 
ity of the band dispersion, which necessarily arises in the 
momentum-space vicinity of each topological charge, as 
follows from the so-called Atiyah-Bott-Shapiro construc¬ 
tion.!^ This property may be viewed as an emergent low- 
energy chiral symmetry, which is characteristic of Weyl 
metals. 

A nearly universal property of metals is the super¬ 
conducting instability, which always exists, at least in 
nonmagnetic metals, at low enough temperature. From 
this viewpoint, the question of superconductivity in Weyl 
metals comes up naturally. Moreover, it is certain in this 
case that the superconductivity will be unconventional 
since, at least at weak coupling, the pairing must occur 
between states in a single nondegenerate conduction or 
valence band. 

As has been shown in Refs. 1571 and 1551 two distinct su¬ 
perconducting states are possible in Weyl metals in prin¬ 
ciple: a BCS (Bardeen-Cooper-Schrieffer) state, in which 
pairing occurs between momentum eigenstates, related to 
each other by inversion symmetry [assuming time rever¬ 
sal (TR) symmetry is broken, but inversion exists]; and 
an FFLO (Fulde-Ferrell-Larkin-Ovchinnikov)-like state, 
with finite-momentum pairs, where the states on the op¬ 
posite sides of each Fermi surface sheet, enclosing indi¬ 
vidual Weyl nodes, are paired. In the latter case, the 
states which are paired are not related to each other by 
any exact symmetry. What makes such a superconduct¬ 
ing state possible in principle, even at weak (but not in¬ 
finitesimal) coupling, is the low-energy chiral symmetry 
of Weyl metals, mentioned above. Moreover, if the inver¬ 
sion symmetry is violated, and as a result, Weyl nodes 
of opposite chirality are shifted away from each other in 
energy to a significant degree,!^ the FFLO state becomes 
the only superconducting state possible, since the exis¬ 
tence of the BCS state relies on inversion symmetry. 

In this paper we show that, in the presence of inversion 
symmetry, the BCS state has a lower energy, at least for 
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short-range momentum-independent pairing interaction 
and in the class of Weyl metals, considered in this pa¬ 
per. This conclusion differ s from that, made in previous 
work on this subjectwhich argued that the FFLO 
state has a lower energy. We will explain the origin of 
this disagreement below. For other recent work on Weyl 
superconductors, see Refs. HOHlll 

The rest of the paper is organized as follows. In Sec¬ 
tion |ll] we introduce the model of a Weyl metal we will 
use, which is based on the TI-NI (normal insulator) mul¬ 
tilayer heterostructure, introduced by one of us before.l^ 
This model has the benefit of being very simple and thus 
amenable to analytic calculations, yet more realistic than 
a generic low-energy model of a Weyl metal would be. In 
particular, it does not suffer from spurious exact chiral 
symmetry of the simplest “relativistic” low energy mod¬ 
els. We also derive in this section the BCS pairing Hamil¬ 
tonian, the form of which plays an important role in our 
analysis. In Sectio n|III| we analyze the BCS Hamiltonian, 
derived in Section |IH using the standard mean-field the¬ 
ory and demonstrate explicitly that, in the presence of 
inversion symmetry the BCS pairing state has a lower 
energy than the FFLO state. In Section |IV] we discuss 
some of the physical properties of the nodal BCS state 
in more detail, in particular the topological properties 
of the mean-field Hamiltonian and the edge states. We 
discuss the differences between our results and previous 
work in Section [V] and conclude by pointing out that, 
when the inversion symmetry (in addition to TR) is vio¬ 
lated, the BCS state is destroyed while the FFLO state 
may survive. 


II. MODEL AND DERIVATION OF THE BCS 
AND FFLO HAMILTONIANS 

As mentioned above, we will start from the model of a 
Weyl metal, based on the TI-NI multilayer heterostruc¬ 
ture, introduced by one of us before.l^^ The model has 
been described extensively in the literature and here we 
will only mention the most essential points, which are 
necessary to understand what follows. The momentum 
space Hamiltonian, describing the multilayer structure, 
is given by 

Hq = vft^{z y. cr) ■'k + i(kz)- (1) 

Here z is the growth direction of the heterostructure, 
(T are Pauli matrices, describing the real spin degree of 
freedom, while t is the pseudospin, describing the top 
(T) and bottom (B) surfaces of the TI layers in the het¬ 
erostructure, and % = 1 units are used here and through¬ 
out the paper. The operator t{kz) describes the motion 
of the electrons in the growth direction of the structure. 
Explicitly it is given by 

i{k.) = tsT^ + Y + h.c) , (2) 

where ts F, are amplitudes for tunnelling between top and 
bottom surfaces of the same (S) or neighbouring (D) TI 


layers and d is the superlattice period. We will take both 
ts,D to be positive for concreteness, this choice does not 
affect any of the physics. This structure describes a Dirac 
semimetal when ts = tu, with two Weyl fermion com¬ 
ponents of opposite chirality residing at the same point 
(0, 0,7r/d) in the first Brillouin zone (BZ). To create a 
Weyl semimetal, we need to separate the Weyl fermions 
of opposite chirality in momentum space. This may be 
accomplished by breaking either TR or inversiorP^ sym¬ 
metries. We choose to break TR and (for now) keep in¬ 
version symmetry intact, as this creates the simplest kind 
of Weyl semimetal state with only two nodes. Breaking of 
TR is accomplished by adding a term b to the Hamil¬ 
tonian Eq. Q 

Hq = vft^{z X cr) • k + t{kz) + ba^. (3) 

Physically this may arise, for example, from magnetized 
transition-metal impurities, introduced into the sample. 
Throughout this paper we will assume that undoped het¬ 
erostructure is almost a Dirac semimetal, i.e. |ts — tnl is 
small. 

We now introduce the simplest and most natural (at 
least for phonon-mediated pairing) kind of pairing inter¬ 
action, i.e. an s-wave short-range, and thus necessarily 
singlet, pairing term, which in second-quantized notation 
has the form 

Hint = -U J 4'|(r)ff'|(r)5'^(r)^'.^(r), (4) 

where U > 0. We want to eventually rewrite Hint in the 
basis of the eigenstates of Hq. To this end, we first write 
the electron field operators in the following way 

vkt (r) = ^ (5) 


Here i labels the unit cells of the TI-NI superlattice in the 
growth direction, a and r are the spin and pseudospin 
labels, (jtiriz) are Wannier-like exponentially-localized 
states, describing the z-axis behavior of the TI surface 
states in the unit cell i and surface t, k = {kx, ky) is the 
transverse momentum, and Lx^y are the sample dimen¬ 
sions in the x- and y-directions. Due to the exponential 
localization of the Wannier functions (j)ir{z)^ it is easy to 
see that, upon substitution of Eq. © into Hint, the dom¬ 
inant terms will correspond to pairing interaction that is 
local in both the i and the t indices 


Hint — 


u 


bi X b / 




EE' 

kk'q ir 


'k+q'it'^ ’ 


( 6 ) 


where 


U = U 



( 7 ) 


Redefining Ud ^ U and transforming to the crystal mo¬ 
mentum basis with respect to the f-index, we then obtain 


H„ 


Y 'y y y y *"k-l-qtT’"k'-gJ,r'"k'ir*"ktTi 

kk'q r 


( 8 ) 
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where k henceforth means the full 3D crystal momentum 
and V = LxLyLz is the sample volume. 

As demonstrated in Refs. 1571 and 1551 the electron pair¬ 
ing will predominantly occur in two distinct channels: 
BCS and FFLO. Correspondingly, we simplify Hint by 
leaving only the two contributions, H^^ and , 

where 


BCS 


TTl 

^int 


^ t t 

kk'r 


( 9 ) 


and 


H 


FFLO 


- y 

V ^ ' 

kk'Qr 


Q+ktr^^Q-k^r^Q-k'ir^^Q+k'tr’ 


( 10 ) 

Here the momentum vector Q labels the locations of the 
Weyl nodes, to be specified below. An important prop¬ 
erty of Eq. (10), which follows from momentum conser¬ 
vation, is that does not couple different Weyl 

nodes. This will play a significant role in selecting the 
lowest energy superconducting state. 

To proceed we diagonalize the noninteracting part of 
the Hamiltonian, Ffo) in order to rewrite the full Hamil¬ 
tonian in the basis of the eigenstates of Hq. Performing a 
canonical (i.e. commutation relations preserving) trans¬ 
formation 


a —)■ r (T 


a T 


( 11 ) 


where The full 4-component eigenvector 

may be viewed as a tensor product 


y{y) = \vr{y)®y{K)). (is) 

The Weyl nodes correspond to points along the z-axis in 
momentum space at which m-(kz) = b — t{kz) vanishes. 
Solving the equation b = t(kz), one obtains Q = Qz, 
where 

Q = .rcc« {yy) 3 5 ± * 0 . (19) 

When the Fermi energy is sufficiently close to the Weyl 
nodes, namely when ep <C b, the Fermi level only crosses 
the s = -I-, r = — band, assuming ep- > 0 for concreteness. 
Projecting onto this band, we write the electron creation 
operators as 

cy = z:t-{yU.^zuyl m 

i.e. we will omit the explicit s = -I-, r = — indices from 
now on in all the equations for brevity. Substituting this 
into and . we finally obtain the following 

projected low energy BCS and FFLO Hamiltonians 

Hbcs = (21) 


which corresponds to the following unitary transforma¬ 
tion of the electron creation operators 

''ktT '^ktTi •'kiT '^kiT> 

•^ktS '^ktS’ ‘^kiB “'^kiBi (^2) 

one obtains 


Ho = vf{z X (t) • k -1- m{kz)a^^ (13) 

where rh{kz) = b -\- t{kz)- The pairing terms are clearly 
unchanged by this transformation. 

Diagonalizing the 2x2 matrix m(kz) one obtains the 
following eigenvalues 


mr{kz) = b + rt{kz), (14) 


where r = ± and t{kz) = \/ts + tjj + 2tstD cos{kzd). 
The corresponding eigenvectors are 


\ugkz)) 


1 ts + tDe-^^^‘^\ 

72 gy )' 


(15) 


Diagonalizing the remaining 2x2 blocks of Hq one finally 
obtains its eigenvalues 


esr(k) = ser(k) = s^v\(kl + kg) + ml{kz), (16) 
where s = ± and the corresponding eigenvectors 


y{y) 


1 

71 



mr(kz) 
egkz) ’ 


'1 - s- 


-{kz) 


gkz) 


( 17 ) 


and 


Hfflo = y^C(k)4ck 

k 

~ ^ /kQ/k'QCQ+kCQ_k'^Q-k''^Q+kO (22) 

kk'Q 


where in Eq. (22) we assume that Qd 1, i.e. Weyl 
nodes are far away from the BZ boundary. 7^) “ ^(J^) ~ 
ep is the band energy, counted from the Fermi energy 
and 


/k 


f m?(kz) 


(23) 


while 


/kQ 


I m{Q + kz) [ m{Q - kz) 
^ [Y e(Q + k) Y e(Q-k) 


+ 



m{Q - kz) m{Q + kz) 
e(Q - k) Y e(Q + k) 


(24) 


The nontrivial momentum-dependent form-factors /k 
and /kQ are a consequence of projection of the pair¬ 
ing interaction terms onto a nondegenerate band. Their 
properties play an important role in the physics of the 
BCS and the FFLO states. It is easy to see that the 
BCS form-factor vanishes identically everywhere on the 
kx = ky = 0 line in momentum space. The four points. 
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at which this line intersects the two Fermi surface sheets 
will produce four nodes in the BCS gap function, as will 
be seen below. On the other hand, the FFLO form-factor 
/kQ never vanishes, which is directly related to the fact 
that the function m{kz) changes sign at the Weyl node 
locations kx = ky = 0, kz — Q. This means that the 
FFLO state is nodeless. It would then seem natural if 
the fully-gapped FFLO state would have a larger conden¬ 
sation energy and thus a lower overall energy than the 
nodal BCS state, with the gap vanishing at four points 
in the BZ. Surprisingly, as we show below, this is not the 
case: the BCS state in fact has a lower energy. 


III. CONDENSATION ENERGIES OF THE BCS 
AND THE FFLO STATES 


It is clear from Eq. (28) that the quasiparticle energy 
indeed vanishes at four points in the first BZ, at which 
the form-factor /j. vanishes. 

The T — 0 BCS equation may be easily solved follow¬ 
ing the well-known steps.l^ One obtains 


Here wd is the Debye frequency, 

= j 

•IT j d 

dkzQ[eF - \m{kz)W, (31) 

-7r/tZ 

is the density of states at Fermi energy, and 


ep’ 




In this section we will evaluate and compare the con¬ 
densation energies of the BCS and the FFLO states, us¬ 
ing the low-energy Hamiltonians, derived in the previous 
section. 


A. BCS state 

We analyze Hbcs using the standard mean-field the¬ 
ory. The mean-field BCS Hamiltonian has the form 


1 f d^k 

is the Fermi surface average of the BCS gap function 

l/kp. 

In the limit ep ^ b, i.e. when the Fermi energy is 
close to the Weyl nodes, the Fermi surface average may 
be easily evaluated explicitly. Indeed, in this case the 
band dispersion in the z-direction in momentum space 
may be assumed to be linear, which follows from the 
leading-order expansion of rn{kz) near the nodes 


Hbcs = ^ 


k 


A 


C(k)CkCk - ^(/kCkclk + /kC_kCk) 


V 


where the pairing amplitude A is given by 


+ ) m(kz)^m(Q) + 

(25) 

where 


dm(kz) 

dkz 


Skz = zLvpSkz, (33) 


k^=Q 


^ ^ H /k(CkC-k) = ^ H /k(C-kCk)- (26) 

k k 

It is important to note that our BCS state has odd parity 
since /k changes sign under inversion. This is already in 
contrast to Ref. EH which claimed an even-parity BCS 
state. 

Diagonalizing Hbcs using Bogoliubov transformation, 
one obtains 

Hbcs = E ^^(k)4V’k + ^ EK(k) " ^(k)] + 

k k 

(27) 

Here "0^ are the Bogoliubov quasiparticle creation oper¬ 
ators, 


VF = Yi,^[b^-its-tDY][{ts + tDY-b% (34) 

is the z-component of the Fermi velocity near the nodes. 
Then one obtains 


9Yf) = 


0 9- 1 

TT^VpVF 


(35) 


and 


(LfkP) = i 

_ I 
“ 4 


1 - 


1 - 


{ni^ikz)) 


vl 

2Cj’ J—ep/vp 


dkz kz 


’ = 6 - <="=> 


The pairing amplitude is thus given by 


E{\Y = x/e(k) + |/kPA2, (28) 


A « (37) 


is the quasiparticle energy and the pairing amplitude A 
satisfies the standard BCS equation (assuming tempera¬ 
ture r = 0) 


To evaluate the condensation energy we take the ex¬ 
pectation value of Hbcs a.t T — 0. In this case there are 
no quasiparticles and we obtain 


U_ I'Ckl^ 

2V ^ E(k)' 

k 


(29) ^=2^EK(k)-^^(k)] + |7- (38) 

k 
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This is again evaluated in the standard waj® and gives 
the following result for the condensation energy, i.e. the 
energy gain compared to the normal state energy Eq 


Ebcs — Eq 
V 


-^5'(ei^)(l/k| 


')A^ 


2bLu'jjg{ep) _i2/[/g(ep) 

" 24 


B. FFLO state 


(39) 


The FFLO state is analyzed in exactly the same way 
as above. One important point to note is that the FFLO 
Hamiltonian Eq. (221 clearly does not mix different Weyl 


nodes, i.e. different values of Q, as required by momen¬ 
tum conservation. This means that, in mean-field theory, 
the contributions of different Weyl nodes may be ana¬ 
lyzed separately and simply summed when calculating 
the total condensation energy. 

The mean-field Hamiltonian for a single Weyl node, i.e. 
with a fixed Q, is given by 


H 


Q 

FFLO 


= E 


^ (/i 


kQ'^Q-l-k'^Q-k 


k 

+ /kQCQ_k‘^Q-rk 

where the pairing amplitude is 

U_ 

V 


V 

W 




(40) 


^ — T A E /kQ(cQ+kCQ-k) 


= ^E'^kQ( 


'^Q-k'^Q+k/- 


Diagonalizing Eq. p0| one obtains 


(41) 


H 


Q 

FFLO 


= E 


k 


e(Q + k)-e(Q-k) 


+ E(k) 




+ i^K(k)-F;(k)] + ^A2, 


(42) 



C(Q + k)+e(Q-k) 


l,/kQPA2. (43) 


A crucial difference between the EFLO and the BCS 
states is the term [e(Q -f k) — e(Q — k)]/2 in the energy 
of the FFLO Bogoliubov quasiparticles in Eq. (42). This 


term is a consequence of the fact that the states with 
momenta Q -I- k and Q — k are not related by any exact 
symmetry, unlike k and —k in the BCS case, which are 
related by inversion. An important property of a Weyl 
metal, however, which may allow the FFLO state to ex¬ 
ist in principle, is the low-energy chiral symmetry, which 
emerges as cf 0 and which implies that 


e(Q + k)«e(Q-k), 


(44) 


the equality becoming more and more precise as the en¬ 
ergy is lowered. Let us first assume that Eq. (44) holds 


exactly, as it would in a low-energy model of a Weyl metal 
with an exactly linear dispersion. Then we obtain 

^FFLO ~ ^ E[^(i^) - ^(k)] + ^A^. 

k k 

(45) 

The mean-field equation for the pairing amplitude, again 
assuming Eq. (44) holds, takes the form, identical to the 
BCS case 


i = —E 

2V ^ 


l./kQp 

£;(k) 


(46) 


To proceed, we now evaluate explicitly the FFLO gap 
function. After elementary algebra, we obtain 


l/i 


kQl 


= F 1 - 


m{Q + kz)m{Q - kz) 


+ A 1- 


'm?{Q + kz) 


' 1 - 


m^{Q - kz) 


Low energy chiral symmetry implies that 
m{Q + kz) « -m{Q - kz)- 
Assuming this to hold, we finally obtain 

l/kQp = 


.(47) 


(48) 


(49) 


which is by a factor of 3/2 larger than the correspond¬ 
ing result for the Fermi surface average of the BCS gap 
function, Eq. (36). This is to be expected, since the 


FELO gap function, unlike the BCS one, does not have 
nodes. This result might then suggest that the EELO 
state should have a lower energy. However, this turns 
out not to be the case. The culprit is another factor 
that influences the magnitude of the pairing amplitude, 
i.e. the density of states. Since in mean-field theory of 
the FFLO state the two Weyl nodes are completely de¬ 
coupled from each other, only half of the total density 
of states determines the magnitude of the gap for each 
Fermi surface sheet. We obtain 

A = (50) 

The extra factor of 2 in the expression above, compared 
to the corresponding Eq. (30) in the BCS case, arises 


precisely from the fact that the density of states per Weyl 
node is g{€F)/2, giep) being the total density of states. 
It is then clear that the magnitude of A in the EELO 
state is exponentially smaller than in the BCS state in 
the weak coupling regime Ug{eF) 1. 

Summing the identical contributions from the two 
Weyl nodes, the total FFLO condensation energy is fi¬ 
nally given by 


= -^5(eF)A2 = -u}l,g{eF)e 


-le/Ugitp) 

(51) 
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Comparing to the corresponding result in the BCS case, 
Eq. (39), it is clear that the FFLO state condensation en¬ 


ergy is exponentially smaller in the weak coupling regime, 
within our model of a Weyl metal. Since the critical tem¬ 
perature Tc ~ A at weak coupling, this also implies that 
Tc of the BCS state is higher than of the FFLO state 
and we thus do not expect any transitions between them 
as a function of temperature. 

In reality the situation is even worse for the FFLO 
state, however, since we have so far been assuming ex¬ 
act chiral symmetry, expressed by Eq. (44). But chiral 


symmetry is not exact, and it is instructive to work out 
the consequences of that. While we have already shown 
that, in our model, the FFLO state is not favored even 
when the chiral symmetry exists, this result is based on 
comparing energies, and thus may not be universal. The 
arguments presented below are of a more general validity. 

Let us expand e(Q±k) in Taylor series with respect to 
the deviation from the Weyl node location k, assuming 
k is small. Let us also take k = for the sake of 
simplicity. We obtain 


e(Q + fc.)-e(Q-fcz)« ^ 


1 (fe 

2 


ki-- 


1 (fe 


kz — Q~\~ 


2 dkl 


kz —Q+ 

kt^ 


de 

dk. 


k^=Q- 

(52) 


k^=Q- 


where we have taken into account the fact that the band 
velocity de/dkz is discontinuous and changes sign at the 
Weyl node locations. Since 


de 

dkz 


kz —Q+ 


de 

dkz 


(53) 


kz=Q- 


the first order term in the expansion above vanishes, 
which is precisely the expression of the approximate chi¬ 
ral symmetry. The quadratic term, however, does not 
vanish and is given by 


1 d^e 

1 dh 

1 

to 

2 dkl 

kz=Q+ ^ 

kz=Q- 


kz=Q 


= d 


2 5 ^- ( t |- t |,)2 


463 


d^h. 


(54) 


We may expect the FFLO state to exist as a local mini¬ 
mum of the free energy as long as 


max |e(Q 

k 


k)-e(Q-k)| < A, 


(55) 


where the maximum is taken over states on the Fermi 
surface. Taking kz ~ ep/vF, where vp dts is the 
Fermi velocity at the Weyl nodes, and assuming b 3> 
I^S — tp], Eq. (55) becomes 




(56) 


where tg should be regarded as the highest energy scale in 
the problem, of the order of the total bandwidth. Using 
Eqs. (35) and (50), we may rewrite this inequality as 


UJD tl ~ Ue% 


(57) 


Using vp ^ vp ^ tsd, it becomes clear that this in¬ 
equality is very hard, if not impossible, to satisfy for any 
reasonable values of the relevant parameters. 

The above arguments lead us to the conclusion that, 
when the inversion symmetry is present and when the 
coupling is weak, nodal odd-parity BCS superconducting 
state, which is a close analog of the A phase in 3He, is 
much more likely to be realized in a Weyl metal. This 
state is topologically nontrivial and is characterized by 
an interesting edge state structure, in which Fermi arcs 
of the “parent” nonsuperconducting Weyl semimetal co¬ 
exist with Majorana edge states of the nodal BCS super¬ 
conductor, as will be described in detail in the following 
section. 


IV. MAJORANA AND FERMI ARC EDGE 
MODES IN THE NODAL BCS STATE 


In this section we will discuss in some detail the non¬ 
trivial momentum-space topology of the nodal BCS state 
and the corresponding Majorana and Fermi arc edge 
states. Some work on this has already been done be¬ 
fore in Refs. sni and SH and the discussion below mostly 
serves to provide a clear connection between these previ¬ 
ously published results and place them in the context of 
our model. 

To discuss momentum-space topology, which involves 
not only states near the Fermi energy, it will be neces¬ 
sary to consider the BCS Hamiltonian without making 
the projection onto the low-energy states only. We will 
still restrict ourselves to the r = — states, since these are 
the bands that touch at the Weyl nodes in the nonsuper¬ 
conducting state, but will take into account both bands 
that touch. The mean-held BCS Hamiltonian then takes 
the form 

H — ^[z;f(z X cr) • k -I- m{kz)<J^ — 

k 

~ 2 ^ ^ ''-kt'^kt) ’ 

k 

where m(kz) = m-{kz), 

^ (59) 

k 


and spin indices were suppressed in the hrst line of 
Eq. ( [^ for brevity. To analyze Eq. (58) we introduce a 
Nambu spinor 


V'k = (Ckt:Ck^>cLkt>cLkt) = (^klt>kU>k2t,'*/'k2t)- 

(60) 

In the Nambu spinor notation, the BCS Hamiltonian 
takes the following form 

77 = i ^[nF(z X cr) • k -I- m(A;^)tT^ - eFK^'jV'kV'k 

k 

- ^ + ^*«”)f^>kV'k- 

k 


(61) 
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FIG. 1. (Color online) Numerically calculated eigenstate dis¬ 
persions for Eq. (611 along the 2 -direction in momentum space 
for a sample of finite width in the ^-direction. Parameters in 
Eq. (611 are chosen as follows: ts = l,to = 0.9, b = l,eF = 
A = 0.2. Majorana modes are pinned to zero energy and exist 
in intervals of width of order |A|/i;F around each Weyl node. 
Fermi arc edge modes are split into particle-hole antisymmet¬ 
ric and symmetric branches at Acf. 


where k is the Nambu pseudospin. 

It is instructive to start from the limit of 6^? = 0, in 


which case Eq. (61) reduces to the model of Ref. flOl In 


this case it is clear that the Nambu pseudospin block of 
the Hamiltonian may be diagonalized separately and we 
obtain 


H = X cr) -k 

k 


-l-p|A|/2]cr^}V'kV'k> 

(62) 

where p = ± labels the two eigenvalues ±|A| of the ma¬ 
trix (Ak+ -h A*k“)/ 2. Apart from the factor of 1/2 in 
front, which expresses the doubling of degrees of free¬ 
dom in the Nambu pseudospin notation, Eq. (62) may 


be viewed as the Hamiltonian of two massive 2D Dirac 
fermions with the masses m(fcz) ± |A|/2 varying as func¬ 
tions of parameter k^. Sign change of the Dirac masses 
signals quantum Hall transitions. In particular, when 
A = 0, m{kz) changes sign at the locations of the two 
Weyl nodes, given by Trjd± ko. The topologically non¬ 
trivial range of k^, where m{kz) > 0, corresponds to the 
range of momenta in which chiral Fermi arc edge states 
exist on any sample surface, not perpendicular to the 
2 -axis.^ Due to doubling of the number of degrees of 
freedom in Eq. (62), these chiral Fermi arc edge states 


come in degenerate particle-hole symmetric and antisym¬ 
metric pairs. Turning on a nonzero A, the Dirac mass 
of the particle-hole symmetric states decreases by |A|/2 
while the mass of the particle-hole antisymmetric states 
increases by the same amount. There is thus an interval 
of momenta near each Weyl node location of width 


Sk, = 


1 


arccos 


t% + ti-ib + \A\/2y 


— arccos 


2tstD 

tl + tl-ib-\A\/2r 

2tstD 


^,(63) 

Vp 




FIG. 2. (Golor online) Numerically calculated eigenstate dis¬ 
persions for Eq. ( |61[ ) along the y-direction in momentum space 
for a sample of finite width in the x-direction, at several 
values of k^. Parameters in Eq. (62 1 are chosen as follows: 
ts = lyto = 0.9,6 = l,eF = 0.5, A = 0.1. (a) kz is outside 
the Fermi surface. Two pairs (corresponding to two surfaces) 
of chiral Fermi arc modes, crossing at Aep at = 0 are 
visible, (b) kz is just inside the Fermi surface. Two chiral 
Majorana modes are visible now, crossing at fey = 0 at zero 
energy, (c) fez is closer to a Weyl node location. Majorana 
modes directly connect with the Fermi arcs. 


in which only the particle-hole antisymmetric states are 
topologically nontrivial (the last equality holds assuming 
|A| ^ b). These momentum intervals give rise to chi¬ 
ral Majorana edge modes, which do not have particle- 
hole symmetric partners, unlike the “ordinary” Fermi 
arc edge states. The Majorana edge states terminate 
at points at which m{kz) = ±|A|/2. These are point 
nodes in the spectrum of the Bogoliubov quasiparticles. 
Each Weyl node in the normal state thus splits into 






























two Bogoliubov-Weyl nodes in the superconducting state, 
which inherit the chirality of the parent Weyl node, as 
seen from Eq. (621. 

When tp > 0, Eq. (61) no longer has the simple form 
of two independent massive Dirac Hamiltonians, but the 
eigenstate spectrum is still easily found. We obtain 


E|(k) = e2(k) + ^ + ± ^4e2(k)e| + |A|2m2(fc,). 

(64) 

The node locations are now given by the solutions of the 
equation 


\m{k^)\ = (65) 

The edge state spectrum may be easily found numeri¬ 
cally, as shown in Eig.[^ The degeneracy of the Fermi arc 
doublet is split, with the particle-hole symmetric branch 
moving down to —ep in energy, while the particle-hole 
antisymmetric branch moves up to ep. The Majorana 
states, which extend between the nodes, remain pinned 
at zero energy (note that zero is at the Eermi energy 
here). 

In the limit |A| ep, we may obtain a simple pic¬ 
ture of the Majorana states from the BCS Hamiltonian, 
projected onto the low-energy s = -I- states 




«k)K“--(A/iK+ + A-/kK-) 




( 66 ) 

where tpu = (cj^, c^k) projected Nambu spinor and 

^(k) = + ky) + m?{kz) — ep. This again has the 

form of the Hamiltonian of a 2D Dirac fermion, with the 
mass ^(0, 0, fcz), which depends on k^ as a parameter. 
The mass changes sign at points, satisfying the equation 
\'>Ti{kz)\ = ep, which coincides with the locations of the 
nodes of the superconducting gap function /k. The topo¬ 
logically nontrivial momentum range, which in this case 
corresponds to negative Dirac mass, coincides with the 
range of momenta inside the two Fermi surface sheets, 
enclosing the Weyl nodes. There are thus Majorana zero- 
energy edge states, which exist on arcs, connecting the 
gap function nodes of the opposite sides of each Fermi 
surface sheet. For sample surfaces perpendicular, say, 
to the ^-direction, the dispersion of the Majorana edge 
modes in the j/-direction is chiral and exist for a given 
kz within the energy interval ±|/kA|. Near the center of 
each Fermi surface sheet, i.e. close to the locations of the 
Weyl nodes along the z axis, the superconducting gap is 
largest in magnitude and the Majorana edge state dis¬ 
persion extends all the way to the Fermi surface, where 
it gets reconnected with the Fermi arc states, as shown 
in Fig.[g 

However, there is a subtlety here, since the projected 
low-energy Hamiltonian Eq. (66) does not actually give 
the correct Majorana edge mode dispersion in the y- 
direction. The reason is that the two Bogoliubov-Weyl 



FIG. 3. Same parameters as in Fig. but values of kz are 
taken near (a) and beyond (b) the Weyl node location. Dis- 
apperance of the Fermi arcs necessitates a “flat band” at 
the transition point. The chirality of the Majorana modes 
changes sign at the transition. The transition itself occurs 
exactly at the Weyl node location in the limit A/ep —>■ 0, but 
is shifted slightly away for any finite A. 


nodes on each Fermi surface sheet have the same chiral¬ 
ity, inherited from the Weyl node, enclosed by the Fermi 
surface sheet, as discussed above. This can not be de¬ 
duced from Eq. (66) and requires analysis of the full un¬ 
projected Hamiltonian Eq. (61). Another way to see this 
is to notice that the Majorana modes must connect with 
the Fermi arcs, which are not low-energy modes and re¬ 
quire the full Hamiltonian Eq. (61) to describe. An in¬ 
teresting consequence of this, first described in Ref. SH 
is a topological transition that happens when one crosses 
the Weyl node locations along the z-axis. The Fermi 
arcs must disappear on approaching Weyl nodes, which 
means that the finite-fcz crossing points in the edge state 
dispersions in Fig.|^must disappear. This leads unavoid¬ 
ably to the appearance of a dispersionless “flat band” at 
the value of kz at which the Fermi arcs disappear, as 
shown in Fig. |3(a)[ As kz is further increased, the Majo¬ 
rana dispersion again becomes chiral, but with chirality 
of opposite sign, as seen in Fig. |3(b)[ 


V. DISCUSSION AND CONCLUSIONS 

We will start this section by discussing differences be¬ 
tween our results and the ones obtained previously in 
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Refs. [571 and [551 which both claimed that the FFLO state 
has a lower energy. In Ref.|38]a specific low-energy linear- 
dispersion model for a Weyl metal was used, in which the 
BCS interaction, Eq. (|^ vanished identically. This is cer¬ 
tainly possible in some realizations of Weyl metals, but 
can not be a general feature. Moreover, such an exact 
cancellation would presumably not happen even within 
the model of Ref. [5H1 if chiral symmetry violating correc¬ 
tions to the low energy model, which are always present, 
as discussed above, were included. However, in agree¬ 
ment with our results, all superconducting states, found 

in Ref. [351 had odd parity. 

Ref. 1371 considered a model, quite similar to ours as far 
as its symmetry properties are concerned, also assuming 
a strictly linear band dispersion in all calculations. How¬ 
ever, the BCS state, found in this reference, was claimed 
to have even parity, which disagrees with our results. The 
FFLO state in this reference also appears to be different, 
with a single order parameter for both nodes, while in our 
case there are two independent (in mean-field theory) or¬ 
der parameters. In fact, a single order parameter would 
be impossible in our case: this would imply internode 
pair scattering already at the level of mean-field theory, 
which violates momentum conservation, as clearly seen 
from Eq. (101. 

We will conclude by pointing out one possible situ¬ 
ation in which the EELO, rather than the nodal BCS 
state, may be realized. FELO state may happen to be 
the ground state if the inversion symmetry, assumed to 
be present in the calculations above, is violated, in ad¬ 
dition to TR symmetry. The presence of at least one 
of those symmetries, i.e. TR or inversion, is necessary 
for the existence of the BCS superconducting state, since 
only those symmetries guarantee that the band eigen¬ 
states at momenta k and — k have the same energy. In a 
Weyl metal, either TR or inversion must be violated, to 
remove the two-fold Kramers degeneracy. In the model 
that we have discussed above, TR was violated from the 
start. This already puts BCS-type superconductivity un¬ 
der some strain, which is manifest in the gap function 
having nodes. When inversion symmetry is violated as 
well, the Weyl nodes will generally be shifted to different 
energies, which implies that the states with momenta k 
and — k will have different energies. Once the inversion 
breaking is strong enough, such that the energy differ¬ 
ence between the Weyl nodes is comparable to the BCS 
pairing amplitude A, given by Eq. (371, the BCS state 
will be destroyed (possibly going through an intermediate 
“ordinary” small-wavevector FFLO state, which we will 


not discuss here). However, the FFLO state is largely un¬ 
affected by this, since its existence relies not on the exact 
microscopic inversion symmetry, but rather on the low- 
energy chiral symmetry, which is unaffected by the bro¬ 
ken inversion symmetry, unless inversion breaking is so 
strong that the energy difference between the Weyl nodes 
becomes comparable to the Fermi energy. Thus, provided 
the inequality Eq. (571 may be satisfied (which, as dis¬ 
cussed above, appears to be hard in the weak coupling 
regime, but may be possible at intermediate coupling), 
the FELO state will be realized. However, it is impor¬ 
tant to keep in mind that since nonlinear chiral symmetry 
violating corrections to the band dispersion are always 
present in a real Weyl metal, realizing the FELO state 
certainly requires a finite pairing interaction strength, as 
the logarithmic divergence of the FFLO pairing suscep¬ 
tibility will always be cut off by the nonlinearity. From 
this viewpoint, a Weyl metal, in which both TR and in¬ 
version are violated, is an example of a metal without 
any weak-coupling BCS instability and thus with a true 
Fermi liquid ground state at T = 0. Another closely re¬ 
lated and interesting possibility, worthy of further, more 
detailed study, is a state, in which the superconducting 
order parameter is nonzero on one of the Fermi surface 
sheets, but is zero on the other one, since the inequality 
Eq. (57) may be satisfied for one of the sheets, but not 
the other, when the inversion sy mmetr y is absent. This 
would be a helical superconductor^^^ coexisting with a 
normal Fermi liquid. Helical superconductors have been 
studied extensively in the general context of supercon¬ 
ductivity in noncentrosymmetric materials,!^ and their 
realization in Weyl metals would be of significant inter¬ 
est. 


Our results may be used as a starting point for further 
studies of superconductivity in Weyl metals. In particu¬ 
lar, since the nodal BCS state is topologically nontrivial, 
as discussed in Section |l^ and in Refs. 00] and 021 the 
question of its electromagnetic response seems to be of 
interest. 
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